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Abstract 

Nascap-2k  is  a  modem  spacecraft  charging  code,  replacing  the  older  codes  NA5CAP/GEO,  NASCAP/LEQ,  POLAR, 
and  DynaPAC.  The  code  builds  on  the  physical  principles,  mathematical  algorithms,  and  user  experience  developed  over 
three  decades  of  spacecraft  charging  research. 


Capabilities  include  surface  charging  in  geosynchronous  and  interplanetary  orbits,  sheath  and  wake  structure  and  current 
collection  in  low-Earth  orbits,  and  auroral  charging.  External  potential  structure  and  particle  trajectories  are  computed 
using  a  finite  element  method  on  a  nested  grid  structure  and  may  be  visualized  within  the  Nascap-2k  interface.  Space 
charge  can  be  treated  either  analytically,  self-consistently  with  particle  trajectories,  or  by  importing  plume  densities  from 
an  external  code  such  as  EPIC  (Electric  Propulsion  Interactions  Code).  Particle-in-cell  (PIC)  capabilities  are  available  to 
study  dynamic  plasma  effects. 

Auxiliary  programs  to  Nascap-2k  include  Object  Toolkit  (for  developing  spacecraft  surface  models)  and  GridT oot  (for 
constructing  nested  grid  structures  around  spacecraft  models). 

The  capabilities  of  the  code  are  illustrated  by  way  of  three  examples:  charging  of  a  geostationary  satellite,  self-consistent 
potentials  for  a  negative  probe  in  a  LEO  spacecraft  wake,  and  potentials  associated  with  thruster  plumes. 


Introduction 


Designers  of  spacecraft  for  government,  commercial,  and  research  purposes  require  advanced  modeling  capabilities  to 
guide  the  design  of  satellites  that  can  survive  and  operate  properly  in  the  natural  environment.  Computer  modeling  of 
flight  experiments  (including  SCATHA,  the  SPEAR3  2  series  and  CHAWS1)  demonstrated  excellent  ability  to  predict 
both  steady-state  and  dynamic  interactions  between  high-voltage  spacecraft  and  the  ambient  plasma.  This  ability  was 
extended  to  inherently  dynamic  problems  involving  three-dimensional  space  charge  sheath  formation,  current  flow  in 
the  quasi-neutral  presheath,  breakdown  phenomena,  plasma  kinetics,  ionization  processes,  and  the  effect  of  unsteady 
processes  on  spacecraft  charging. 

NASCAP/GEO4  5  6  (NASA  Charging  Analyzer  Program  for  GEosynchronous  Orbit)  was  the  standard  tool  for  the 
computation  of  spacecraft  charging  in  tenuous  plasmas  for  more  than  two  decades.  Since  then,  the  fully  three- 
dimensional  computer  codes  NASCAP/LEO7,8  (NASA  Charging  Analyzer  Program  for  Low-Earth  Orbit),  POLAR9 
(Potentials  Of  Large  objects  in  the  Auroral  Region),  and  DynaPAC 30  (Dynamic  Plasma  Analysis  Code)  were  developed 
to  address  various  other  spacecraft-plasma  interactions  issues.  Nascap-2kn  builds  on  the  capabilities  of  these  older 
codes,  giving  the  spacecraft  designer  much-improved  modeling  capabilities  by  taking  advantage  of  a  greater 
understanding  of  the  pertinent  phenomena,  employing  more  advanced  algorithms,  and  implementing  a  state-of-the-art 
user  Interface,  including  three-dimensional  post-processing  graphics.  Nascap-2k  is  being  developed  as  part  of  a  program 
sponsored  jointly  by  the  Air  Force  Research  Laboratory  at  Hanscom  AFB  and  by  NASA's  Space  Environments  and 
Effects  (SEE)  Program  at  Marshall  Space  Flight  Center.  The  current  release  is  Version  3.0. 

Nascap-2k  is  an  interactive  toolkit  for  studying  plasma  interactions  with  realistic  spacecraft  in  three  dimensions.  As  it 
incorporates  physics  developed  for  all  the  previous  codes,  it  can  solve  problems  appropriate  to  both  tenuous  (e  g.,  GEO 
orbit  or  interplanetary  missions)  and  dense  (e.g.,  LEO  orbit)  plasma  environments.  Nascap-2k  is  targeted  to  spacecraft 
design  engineers,  spacecraft  charging  researchers,  and  aerospace  engineering  students.  The  graphical  user  interface  is 
designed  to  help  less  experienced  users  easily  solve  moderately  complex  plasma  interactions  problems.  Figure  1  shows 
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several  views  of  the  Nascap-2k  interface,  including  the  main  problem  setup  page,  the  GridTool  interface,  and  views  of 

particle  trajectories  and  space  and  surface  potentials. 

The  core  capabilities  of  Nascap-2k  include: 

1 ,  Define  spacecraft  surfaces  and  geometry  and  the  structure  of  the  computational  space  surrounding  the  spacecraft; 

2,  Solve  for  time-dependent  potentials  on  spacecraft  surfaces; 

3,  Solve  the  electrostatic  potential  about  the  object,  with  flexible  boundary  conditions  on  the  object  and  with  space- 
charge  computed  either  fully  by  particles,  fully  analytically,  or  in  a  hybrid  manner; 

4,  Generate,  track,  and  otherwise  process  electrons  and  ions,  represented  as  macroparticles  in  the  computational  space; 
and 

5,  View  surface  potentials,  space  potentials,  particle  trajectories,  and  time-dependent  potentials  and  currents. 

To  accomplish  these  capabilities,  Nascap-2k  consists  of: 

l  The  Nascap-2k  graphical  user  interface  (Figure  1),  a  user-friendly  environment  for  definition  of  problem 
parameters,  strategizing  and  running  calculations,  and  visualizing  results; 

2.  Object  Tool  Kit  (OTk),  an  interactive  program  for  the  definition  of  spacecraft  surfaces; 

3.  GridTool ,  an  interactive  program  to  create  arbitrarily  subdivided  grids  in  the  space  surrounding  a  spacecraft  model; 

4.  The  modules  that  comprised  the  DynaPAC  code  as  the  major  computational  engine.  These  have  been  converted  to 
DLLs  (dynamic  link  libraries)  to  run  seamlessly  within  Nascap-2h ; 

5.  A  new  analysis  module  implementing  the  Boundary  Element  Method12  (BEM)  for  calculating  surface  charging  in 
Geosynchronous  Earth  Orbit  (GEO),  in  the  Solar  Wind,  or  in  other  tenuous  plasma  environments. 


Figure  I,  Views  of  the  Nascap-2k  tabbed  user  interface,  including  main  problem  setup  (center  left),  particle  trajectories 
(upper  right),  surface  and  space  potentials  (center  right),  and  GridTool  interface  (bottom  center). 
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Nascap-2k  uses  a  high-order  finite  element  representation  for  the  electrostatic  potential  that  assures  that  electric  fields 
are  strictly  continuous  throughout  space.  The  electrostatic  potential  solver  (originally  developed  for  DynaPAC  n)  uses  a 
conjugate  gradient  technique  to  solve  for  the  potentials  and  fields  on  the  spacecraft  surface  and  through  the  surrounding 
space*  Several  analytic  and  numerical  space  charge  density  models  are  available,  including  Laplacian,  Linear,  Non¬ 
linear,  Frozen  Ions,  Full  Trajectory  Ions,  Hybrid  PIC  (appropriate  to  the  several  microsecond  timescale  response  to  a 
negative  pulse),  and  Full  PIC. 

Particle  tracking  is  used  to  study  sheath  currents,  to  study  detector  response,  to  generate  steady-state  charge  densities,  or 
to  generate  space  charge  evolution  for  dynamic  calculations.  Naseap-2k  generates  macroparticles  (each  of  which 
represents  a  collection  of  particles)  either  at  a  “sheath  boundary”,  the  problem  boundary,  or  throughout  all  space. 
Alternatively,  particles  can  be  initialized  with  a  user-generated  file.  Particles  are  tracked  for  a  specified  amount  of  time, 
with  the  timestep  automatically  subdivided  at  each  step  of  each  particle  to  maintain  accuracy.  The  current  to  each 
surface  of  the  spacecraft  is  recorded  for  further  processing. 

The  Results  3D  tab  of  the  Nascap-2k  user  interface  is  used  to  generate  graphical  output  illustrating  such  quantities  as 
object  surface  potentials,  space  potentials,  particle  positions,  or  particle  trajectories.  Using  Java3D  capabilities,  these 
figures  can  be  rotated,  panned,  zoomed,  or  measured.  Contour  levels  and  other  plotting  attributes  are  modi  fed  through 
the  user  interface.  The  Results  tab  is  used  to  view  time  histories  and  obtain  numerical  values  for  potentials  and  surface 
currents. 

The  modular  structure  of  Nascap-2k  is  illustrated  in  Figure  2*  Surface  charging  is  done  in  the  new  BEM  module.  Space 
potentials  and  particle  trajectories  are  calculated  with  DLLs  built  from  the  DynaPAC  modules10  34  The  suite  of  codes  is 
written  in  Java  (user  interface),  C++  and  Fortran  (science),  and  C  (utility  routines)  and  is  maintained  on  the  Win32 
platform. 


Figure  2,  Nascap-2k  Module  Diagram. 
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Spacecraft  Models  and  Object  Toolkit 


A  Nascap-2k  application  usually  begins  by  building  a  geometrical  spacecraft  model  with  Object  Toolkit.  Objects  arc 
built  using  the  five  native  components  shown  in  Figure  3,  together  with  components  imported  from  standard  finite 
element  preprocessing  software  and  components  previously  defined  and  saved  using  Object  Toolkit.  Direct  mesh  editing 
capabilities  are  available  to  create  subdivision  and  to  build  components  with  complex  shapes.  Figures  4  through  6  show 
some  examples  of  spacecraft  defined  using  Object  Toolkit . 


Figure  3,  Native  components  (clockwise  from  upper  left:  Box,  Boom,  Panel,  Cylinder,  Dish)  available  within  Object 
Toolkit  for  defining  spacecraft  models. 

Bach  elemental  surface  of  a  spacecraft  model  has  attributes  of  “Conductor  Number"  and  “Material  Name."  The 
conductor  number  attribute  is  used  to  represent  electrical  circuitry  coupling  the  spacecraft  surfaces,  such  as  biasing  of 
surfaces  or  capacitive/resistive  coupling.  With  each  material  name  is  associated  a  list  of  properties  including  thickness, 
bulk  and  surface  conductivity,  and  photoemission  and  secondary  electron  emission  coefficients.  Material  properties  can 
be  edited  in  both  Object  Toolkit  and  the  main  Nascap-2k  interface.  In  addition,  advanced  features  such  as  grounding  tabs 
can  be  placed  on  the  spacecraft  model. 
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Figure  4.  Simple  spacecraft  model,  featuring  a  solar  array  yoke  (upper  right)  built  in  Object  Toolkit  from  three  booms. 


Figure  5.  Nascap~2k  model  of  the  MESSENGER  spacecraft,  showing  biased  solar  array  surfaces. 
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Figure  6.  DMSP  spacecraft  model  constructed  in  Object  Toolkit. 
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Nascap-2k  Papers  and  Examples  at  SCTC  Conferences 


Descriptions  of  Nascap-2k  and  its  predecessor  codes,  and  simulations  and  analyses  based  on  those  codes,  have  been 
presented  at  every  Spacecraft  Charging  Technology  Conference  dating  back  to  the  first  conference  in  Colorado  Springs, 
October  1976.  Table  1  shows  a  list  of  the  papers,  together  with  the  examples  contained  in  each,  for  the  three  most  recent 
conferences. 


Table  l.  Nascap-lk  papers  presented  at  re  cent  Spacecraft  Charging  Technology  Conferences. 


9th  Spacecraft  Charging  Technology  Conference  (Tsukuba,  Japan,  April  2005) 

■ 

Nascap-2k  Spacecraft  Charging  Code  Overview  (This  Paper) 
o  Geosynchronous  Charging 
o  Current  Collection  in  Wake 
o  Potentials  in  Ion  Engine  Plume 

■ 

Nascap~2k  Simulations  of  a  VLF  Plasma  Antenna 
o  PIC  Capabilities 

• 

Reverse  trajectory  approach  to  computing  ionospheric  currents  to  the  Special  Sensor  Ultraviolet  Limb  Imager 
on  DMSP 

O  Determine  current  to  sensor 

Sth  SCTC  (Huntsville,  Alabama,  October  2003) 

• 

Validation  of  Nascap-2k  Spacecraft-Environment  Interactions  Calculations 
o  Current  collection  by  sphere 

o  Geosynchronous  charging  (comparison  with  other  codes) 

□  Bipolar  sheath  collection  (comparison  with  older  results) 
o  Auroral  charging  (comparison  with  older  results) 

* 

Nascap-2k  -  An  Overview 

o  GEO,  interplanetary,  and  auroral  charging 
o  LEO  charging  showing  vxB 

■ 

Nascap-2k  as  a  PIC  Code 

o  Dynamic  antenna  sheath  calculations 

T"1  SCTC  (Noordwijk,  The  Netherlands,  April  2001) 

■ 

MESSENGER  Spacecraft  Charging  Analysis 

o  Charging  in  Hermean  and  near-Hermean  environment 

* 

Comparison  of  the  NASCAP/GEO,  SEE  Interactive  Charging  Handbook,  and  Nascap-2k  Spacecraft 
Charging  Codes 

o  Geosynchronous  orbit  charging 

■ 

Nascap-2k  Spacecraft  Charging  Models:  Algorithms  and  Applications 
o  Implementation  Details 
o  Geosynchronous  orbit  charging 
o  Solar  Wind  charging 

Page  7  of  21 


Example  I:  Geosynchronous  Orbit  Charging 

The  earliest  and  most  common  application  for  this  family  of  codes  is  to  study  charging  of  spacecraft  in  geostationary 
orbit.  Nascap-2k  is  designed  to  make  this  type  of  analysis  particularly  easy, 

A  simple  Object  ToolkihbuUt  model  of  a  geosynchronous  spacecraft  is  shown  in  Figure  7.  Apart  from  the  solar  arrays, 
the  spacecraft  is  covered  with  low-secondary-emission  insulators,  except  for  some  patches  of  high-secondaiy-emission 
OSRs  on  the  side  of  the  box. 

On  the  initial  Problem  tab  of  Nascap-2k  {Figure  8)  we  load  the  object  and  select  a  surface  charging  problem  in  a 
geostationary  environment.  Since  Nascap-2k  can  calculate  surface  potentials  and  electric  fields  using  the  Boundary 
Element  Method  (BEM),  there  is  no  need  to  define  a  spatial  grid,  and  consequently  all  options  requiring  a  grid  are 
disabled.  With  a  grid,  we  would  be  able  to  calculate  and  display  space  potentials  and  particle  trajectories. 

Next  we  define  the  environment  on  the  Environment  tab,  as  shown  in  Figure  9.  Any  single  or  double  maxwell ian 
environment  can  be  defined,  and  some  predefined  environments  are  provided.  In  this  case,  we  use  the  built-in  “NASA 
Worst  Case”15  environment.  The  magnetic  field  and  the  sun  direction  and  intensity  (1.0  for  sunlit  and  0,0  for  eclipse)  are 
also  set  here.  Note  that  the  tab  as  shown  is  customized  for  geosynchronous  substorm  parameters;  if  a  low- Earth -orb it  or 
auroral  environment  is  selected  the  tab  looks  quite  different. 


mm 


ChefnGlazs 


Figure  7.  Simple  geosynchronous  spacecraft  used  for  example. 
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Figures  10  and  1 1  show,  respectively,  the  Charging  and  Script  tabs  that  customize  how  the  calculation  is  performed.  In 
this  case,  the  timestepping  parameters  were  edited  to  produce  more,  shorter  time  steps  than  the  default,  and  the 
automatically  generated  script  was  used  without  modification.  Exposing  the  script  provides  a  great  deal  of  flexibility  in 
running  problems  with  Nascap-2L  Customized  scripts  can  be  built  and  edited  both  within  the  Script  tab  and  externally 
using  a  text  editor  or  specialized  XML  editor. 

The  calculation  is  initiated  by  clicking  the  “Run  Script"  button.  Monitors,  such  as  the  “BEM  Monitor”  shown  in  Figure  12, 
appear  to  indicate  the  progress  of  the  calculation.  In  this  case,  the  minimum  and  maximum  potentials,  simulation  time,  and 
elapsed  time  are  periodically  updated.  Most  of  the  two  minutes  of  elapsed  time  shown  on  the  monitor  is  attributable  to  the 
calculation  of  the  BEM  matrix  elements,  which  are  stored  on  disk  so  they  need  not  be  recalculated  for  additional 
simulations. 


-  NascapZk  -  E:\NASCAP2000^SCTC2005\N2kOvervTew\Geosynchronous\DSCS\0SCSProject.xml 


Erie  gm  View  Materials  Help 

Probtem  [^Enytronmenl  Applied  Potentials  *  Charging  Space  r^nntials  ]  Hart  teles  |  Script  ■  Results  [  Results  3D 


EE® 


n  Orbit  ■ 


GEO 

Oleo 

O  Polar 

O  Interplanetary 

C  Radiation  Belt 


Object  Status:  Loaded 


Edit  Object.- 


Grid  Status:  NO f  leaded 


Edit  Grid.- 


-  Problem  Type - 

0  Surface  Charging 
(#■  Analytic  Currents 

O  T  lacked  Part  (do  Currents  (Preliminary  model) 
O  Analytic  Electron  &  Tracked  Jon  Oji rents 
O  Custom  Current  Module 
□Potentials  in  Space  nr  Detector  Analysis 
Analytic  Spore  Charge 
O  Self  consistent  utfrUi  ton  trajectories 
O  Consistent  with  EPIC  ton  Densities 
□Surface  Currents 

□  Time  Dependent  Plasma 

□  Deep  Dielectric 


Figure  8.  Opening  Problem  tab  of  Nascop~2k  interface  set  for  geosynchronous  charging  example. 
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Figure  9,  Naseap-2k  Environment  tab,  customized  for  GEO  environment  definition* 


Nascap2k  -  E;\NASCAP2000\SCTC2005\N2kOverview\Geosynchronoys\DSCSVDSCSProject,)tmt 


pe  Edit  view  Materials  tfeip 

Problem  ]  Environment”]  Applied  Potentials  1  Charging  v  Space  Potent ials~|  Particles  |  Scrip!  [  Results  ]  Results  30^ 


-  Charging  time 


Start  Time  (sec):|Q  o| 
Min  Time  step  (sec):!o  1 00 
Number  of  Timesteps*:|90 


End  Time  (sec):|300  0_ 


Max  Tiiuestep  (see):' 10.00 


only  for  Aflstytto  Currents)- 


Figure  10.  Nascap-2k  Oh  urging  tab.  used  to  set  timestep  parameters  for  charging  calculations. 
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Figure  11.  Nascap~2k  Script  tab,  set  for  a  simple  geosynchronous  environment  charging  calculation.  While  this  script 
was  generated  automatically,  more  complex  and  customized  scripts  can  be  built  and  edited  either  within  this  tab  or 
externally  using  a  text  or  XML  editor. 
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Figure  12,  "'BEM  Monitor”  indicates  the  progress  of  a  calculation  using  the  BEM  module,  such  as  a  surface  charging 
calculation. 
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Calculation  results  are  post- processed  and  displayed  on  the  Results3D  tab  (Figure  13)  and  on  the  Results  tab  (Figure  14). 
Figure  13  shows  that,  after  five  minutes  of  charging,  most  of  the  dark  surface  material  charges  to  about  10  kV  negative, 
with  the  emissive  OSR  surfaces  about  8  kV  negative.  The  least  negative  surfaces,  at  about  7  kV  negative,  are  at  the  center 
of  the  sunlit  side  of  the  central  box  and  toward  the  outboard  end  of  the  solar  arrays.  As  shown  in  Figure  13,  we  can  select 
individual  surfaces  to  obtain  detailed  information  about  them.  In  a  gridded  calculation,  this  tab  can  define  and  display 
contour  planes  of  external  potential  and  particle  trajectories  superposed  on  the  surface  plot. 

The  Results  lab  (Figure  14)  is  used  to  display  detailed  numerical  and  time-dependent  information.  The  left  part  of  the 
tab  allows  specification  of  groups  of  surface  cells,  individual  surfaces,  or  conductor  entities,  and  selection  of  potentials, 
electric  fields,  or  currents.  Shown  in  Figure  14  are  the  conductor  potential  (most  negative  curve)  and  the  minimum, 
maximum,  and  average  potentials  of  the  solar  cells.  The  plotted  information  may  be  displayed  as  text  on  the  Text  subtab 
and  pasted  into  a  spreadsheet  for  further  processing.  Figure  15  shows  the  same  data  plotted  in  Microsoft  Excel,  along 
with  the  potential  of  the  dark  insulator  that  covers  most  of  the  satellite  surface. 


Figure  13.  Results3D  tab  showing  surface  potentials  at  end  of  calculation,  and  detailed  information  for  a  selected 
surface  cell. 
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Figure  14,  Results  tab,  showing  time  development  of  ground  potential  and  range  of  potentials  on  solar  array  surfaces. 
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Figure  15.  Charging  results  imported  into  Microsoft  Excel  from  Text  subtab. 
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Example  2;  Self-Consistent  Potentials  in  a  Wake 

To  study  ion  collection  by  a  high  voltage  object  in  a  spacecraft  wake,  the  CHAWS16  (Charging  Hazards  and  Wake 
Studies)  experiment  was  flown  about  the  WSF  (Wake  Shield  Facility)  on  STS-6G  (February  1994)  and  STS-69 
(September  1995).  Originally  simulated  using  the  DynaPAC  code,  we  use  this  problem  as  the  primary'  example  of  a  self- 
consistent  (between  potentials  and  trajectories)  calculation  in  Nascap-Ik. 

Figure  16  shows  the  Object  Toolkit  model  of  the  WSF  with  CHAWS.  The  CHAWS  probe  is  the  off-center  we  1 1 -resolved 
rod  seen  toward  the  left  of  the  figure,  which  shows  the  wake  side  of  the  WSF,  Negative  biases  of  up  to  2  kV  were 
applied  to  the  probe,  which  was  instrumented  to  measure  both  the  total  current  and  the  distribution  of  current  over  the 
surface* 

Since  we  need  to  calculate  both  space  potentials  and  particle  trajectories,  a  grid  structure  is  built  in  the  space 
surrounding  the  model.  The  GridTool  representation  of  the  grid  structure  is  shown  in  Figures  17  and  18,  The  grid  is 
designed  to  provide  high  resolution  both  in  the  neighborhood  of  the  probe  and  in  region  of  plasma  flow  that  is  perturbed 
by  the  applied  potential. 


— [tf  Spacecraft  - 

Display: 

Material 

▼ 

UlcUa** 

view: 

Stanford 

Figure  16*  Geometric  model  of  the  WSF  wake  side,  showing  the  CHAWS  probe  (well-resolved  cylinder  at  left). 
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Figure  17.  GridTool  view  of  grid  structure  used  for  the  CHAWS  problem.  The  right  pane  details  the  grid  and  subgrid  structure. 


Figure  18.  GridTool  view  of  highly  resolved  grid  near  CHAWS  probe. 
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Figure  19  shows  the  environment  specification  for  this  problem.  For  a  spacecraft  orbiting  in  a  LEO  plasma,  we  must 
specify  not  only  the  plasma  density  and  temperature,  but  also  the  spacecraft  velocity  and  the  ion  species,  which  together 
determine  the  "ram  ion  energy.  Figure  20  shows  the  script  for  this  calculation,  which  is  generated  automatically  by 
Nascap-2k.  The  script  directs  the  code  to  generate  and  track  particles  in  the  existing  potentials,  and  then  recalculate  the 
potentials  using  space  charge  derived  from  the  ion  trajectories.  (The  electron  density  is  assumed  to  be  barometric,)  The 
sequence  repeats  to  approach  a  self-consistent  solution.  The  script  also  refers  to  input  files  (which  are  generated 
automatically  and  can  also  be  user-edited)  to  the  Nascap-2k  modules  and  to  the  output  files  from  those  modules. 

Figure  21  shows  a  subset  of  ion  trajectories  for  this  problem.  Ions  are  generated  at  the  problem  boundary1  with  a  thermal 
spread  about  the  ram  direction  and  are  tracked  until  they  strike  the  object  or  leave  the  problem  space.  Ion  space  charge  is 
accumulated  in  the  grid  at  each  tracking  step,  and  the  current  corresponding  to  ions  reaching  to  CHAWS  probe  is 
recorded.  It  is  noteworthy  that  no  ions  strike  the  side  of  the  probe  nearest  the  edge  of  the  WSF,  but  rather  ions  either 
strike  the  tip  of  the  probe  or  else  miss  the  tip  and  are  attracted  back  to  the  side  of  the  probe  facing  the  center  of  the  WSF. 

Finally,  Figure  22  displays  the  self-consistent  potentials  about  the  CHAWS  probe  behind  the  WSF.  Note  that  the 
potential  field  of  the  probe  extends  considerably  into  the  ram  flow  in  the  vicinity  of  the  near  edge  of  the  WSF,  and  this 
influence  manifests  itself  in  the  trajectory  deflection  for  this  region  seen  in  Figure  21.  The  grounded  rear  surface  of  the 
WSF,  the  WSF  instrumentation,  and  the  ion  space  charge  all  play  a  role  in  limiting  the  penetration  of  CHAWS  potentials 
into  the  wake,  so  that  the  potential  Is  shielded  much  faster  than  would  be  the  case  if  the  wake  were  a  cylindrical  vacuum 
region. 
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Figure  19,  Environment  tab,  showing  specification  of  plasma,  spacecraft  velocity,  and  particle  species  for  the  CHAWS 
problem. 
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Figure  20.  CHAWS  problem  script  to  iteratively  create  and  track  ions  and  solve  potentials  in  the  space  charge  of  those 
ions. 
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Figure  21  Trajectories  of  ram  particles,  including  thermal  spread,  for  CHAWS  problem. 
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Figure  22,  Self-consistent  potentials  about  the  CHAWS  probe  in  the  wake  of  the  WSF. 

Example  3:  Potentials  in  Electrostatic  Thruster  Plumes 

A  recent  enhancement  made  to  Nascap-2k  is  the  inclusion  of  ion  densities  from  plasma  sources  (such  as  thruster  plumes) 
in  the  computation  of  potentials  in  space.  One  reason  for  needing  potentials  due  to  plasma  sources  is  the  computation  of 
contaminant  trajectories.  We  plan  to  extend  this  capability  to  study,  among  other  phenomena,  the  influence  of  spacecraft 
surfaces  on  engine  plumes  and  interactions  between  plumes. 

For  this  example  the  ion  density  due  to  plasma  sources  is  imported  from  the  EPIC  (Electric  Propulsion  Interactions 
Code)  computer  code17  which,  like  Noscap-2k>  has  been  developed  by  SA1C  for  NASA's  SEE  program.  The  import  is 
done  via  SOAP  (an  XML-based  interprocess  communication  scheme  formerly  known  as  Simple  Object  Access 
Protocol),  Nascap-2k  then  solves  a  nonlinear  Poisson  equation  for  space  potentials  using  the  surface  potentials  as 
boundary  conditions.  As  was  done  in  the  CHAWS  example  presented  above,  the  charge  density  consists  of  a  known  ion 
density  plus  a  barometric  electron  density  ; 


—  =  —  (l  -expfO-ifrt,  )/©)) 

K=61n(^] 

where  0  and  n  are  the  plume  electron  temperature  and  reference  density  respectively,  p  is  the  total  charge  density  to  be 
used  in  Poisson's  equation,  and  piq„  is  the  known  ion  density,  in  this  case  imported  from  EPIC . 

We  treat  a  configuration  that  might  characterize  a  geosynchronous  spacecraft  during  North-South  station-keeping.  On 
the  same  spacecraft  used  above  for  the  geostationary  charging  example,  we  affix  two  NSTAR18  thrusters  pointing 
Northwest  and  NorthEast,  (A  thruster  as  large  as  NSTAR  is  unlikely  to  be  used  for  this  purpose,  but  w^e  selected  it 
because  we  had  a  good  model  for  its  plume.)  Figure  23  shows  the  spacecraft  and  plasma  densities  as  displayed  in  EPIC. 
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Figure  23.  Plasma  density  displayed  in  EPIC  for  two  NSTAR  thrusters,  as  described  in  text 

For  a  reference  density  we  chose  a  value  of  I  *  1012  m“\  which  is  the  density  of  charge  exchange  ions  just  outside  the 
main  beam.  This  density  is  also  used  as  the  minimum  ion  density,  so  that  regions  void  of  plume  ions  are  set  to  zero 
potential.  The  beam  electron  temperature  was  taken  as  I  eV.  We  set  the  average  potential  on  the  solar  arrays  to  be  +20 
V,  and  the  potential  of  all  other  surfaces  to  be  -2  V,  Figure  24  shows  the  resulting  potentials  in  space. 


Figure  24.  A fa$cap-2k  display  of  space  and  surface  potentials  in  the  presence  of  the  thruster  plumes  shown  in  Figure  23. 
Other  parameters  as  described  in  the  text. 
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Summary 


Nascap-2k  is  a  user-friendly  code  for  the  study  and  analysis  of  spacecraft-plasma  interactions  in  a  variety  of  important 
space  environments.  The  code  builds  on  many  years  of  study  and  experience  in  building  codes  and  performing  analyses. 
Many  examples  of  Nascap-2k  calculations  have  been  presented  at  this  series  of  conferences.  Additional  examples  of 
calculations  by  Nascap-Ik  and  its  predecessor  codes  may  be  found  in  the  published  literature  dating  back  to  the  late 
1970’s.  In  this  paper  we  show  a  quick  walk-through  of  a  geosynchronous  charging  calculation,  and  results  for  potentials 
in  a  spacecraft  wake,  and  for  potentials  in  thruster  plumes. 

Nascap-2k  is  distributed  by  the  Space  Environments  and  E fleets  Program  at  NASA’s  Marshall  Spaceflight  Center, 
Huntsville,  Alabama,  http://see.msfc.nasa.gov/. 
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